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Abstract 

Chemical reactions can be modelled via diffusion processes conditioned 
to make a transition between specified molecular configurations representing 
the state of the system before and after the chemical reaction. In particular 
the model of Brownian dynamics - gradient flow subject to additive noise - 
is frequently used. If the chemical reaction is specified to take place on a 
given time interval, then the most Ukely path taken by the system is a mini- 
mizer of the Onsager-Machlup functional. The F— limit of this functional is 
determined in the case where the temperature is small and the transition time 
scales as the inverse temperature. 



1 Introduction 

In this paper we study the problem of determining the most likely paths taken by 
a system undergoing a chemical reaction. We employ the model known as Brow- 
nian dynamics [IJ: the atomic positions are assumed to be governed by a gradient 
flow in an interaction potential, subject to additive thermal noise. The resulting 
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stochastic differential equation is conditioned to make a transition between two 
different atomic confirgurations representing the state of the system before and af- 
ter the chemical reaction [2] . If the chemical reaction is specified to take place on 
a given time interval, then the most likely path taken by the system is a mmimizer 
of the Onsager-Machlup functional [3, 4|. In |5| numerical computations are pre- 
sented which study minimizers of this functional for a variety of model problems in 
low dimension, together with some higher dimensional problems from physics and 
chemistry such as vacancy diffusion and the Lennard- Jones cluster of 38 atoms. 
The minimizers exhibit a number of interesting effects, including multiple path- 
ways between configurations, together with transition paths which concentrate at 
saddle points of the potential, and not at minima. The computations in [5] were per- 
formed at low temperature over fixed long intervals. This paper is concerned with 
determining the F— limit (see |6, 7|) of the Onsager-Machlup functional in the case 
where the temperature is small and the transition time scales as the inverse temper- 
ature. Minimizers of the Onsager-Machlup functional have received considerable 
attention in the chemistry literature (see tSJ and [5] and the references therein). 
One of the motivations for the work in Q is to show that there are certain artefacts 
in this minimization procedure which could be construed as unphysical. In this 
paper we provide a mathematical theory to explain the computational observations 
in The F— limit which we exhibit may also be of independent interest. The 
F— limit can be optimized via knowledge of the critical points of the potential, the 
trace of the Hessian of the potential at the critical points, and certain heteroclinic 
orbits connecting them. Notably these heteroclinic orbits are in a Hamiltonian sys- 
tem; heteroclinic orbits in the forward or backward gradient flow found from our 
model at zero temperature are solutions of this Hamiltonian problem but not nec- 
essarily vice versa. Regarding the role of the Hessian in the F— limit, it is pertinent 
to mention the paper [9| in which second derivative of the potential plays a role in 
a large deviations principle for SDEs. 

In section |2] we provide a precise mathematical description of the conditioned 
SDE which forms our mathematical model. We provide an informal derivation of 
the Onsager-Machlup functional whose minimizers determine most likely transi- 
tion paths and we explain the sense in which this informal argument can be made 
rigorous. In section[3|we compute the F— limit of the Onsager-Machlup functional 
building on related analyses in f TOl [TTll . The paper concludes, in section |4} with 
an informal characteriztion of the F— limit, together with numerical experiments 
which illustrate this characterization. 



2 Set-Up 

Consider the following conditioned SDE for x £ C([0, T]; M'^) making a transi- 
tion between two states and 2;+ in time T: 

dx = -DV{x)dt + \/2edW, 

_ ' (2.1) 

x(0) = X and x{T) = x . 

Here F : R is the potential, W is a standard Brownian motion in M 

and e <^ 1. In many applications in physics/chemistry N = nd, d G {1,2,3} 
being the physical dimension and n the number of particles. Then e is the (non- 
dimensional) temperature and x G C([0, T]; M^) denotes the configurational path 
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of the n atoms making a transition between . Throughout we assume that 
are chosen as critical points of V. We will be particula rly i nterested in choosing 
these critical points to be minima, so that the problem ( |2.1| describes a chemical 
reaction; but, m the numerical computations, we will also choose saddle points in 
order to illustrate certain mathematical phenomena. 

Escape from local minima occurs on timescales which are exponentially long 

in e~^, whilst transitions themselves occur on scales which are logarithmic in e^^. 
We focus interest on an intermediate timescale between these two regimes, which 

is e~^. This timescale is hence long enough to capture a single transition, but not 
long enough to capture the typical waiting time in a potential minimum. In partic- 
ular, transition paths calculated in this scaling which pass through an intermediate 
minimum x^ of V will not necessarily exhibit the "typical" behaviour that would 

be exhibited if the time to make the transition from x^ to x^ via x^ were left free. 
The behaviour of the most likely single transitions in this regime is captured by 
minimizers of the Freidlin-Wentzell action |'T2irT3'|. Numerical methods to capture 
multiple such transitions are the subject of active study, and the reader should con- 
sult [14J, and the references therein, for details. From an applied perspective, our 
work serves to highlight the potential pitfalls of using the Onsager Machlup ap- 
proach to compute transition paths, and this point is discussed in detail in 15]. As 
mentioned above, the work in this paper gives a mathematical explanation for the 
numerical results observed in [5 |. The results may also be of independent interest 
from the point of view of the calculus of variations. In particular the results show 
how hetereoclinic orbits in a certain Hamiltonian flow form the building block for 
construction of the F— limit. This fact may also be of interest in applications where 
attention has focussed on heteroclinic orbits in the forward or backward gradient 



flow found from ( 2.1[ ) with e = 0; these form particular instances of heteroclinic 
orbits in the Hamiltonian flow, but non-gradient connections are also possible, as 
we will demonstrate. 

To enforce the scaling of interest we choose T = e^^ and rescale time as 
t = e^^s, to obtain 

dx = --DV(x)ds + V2dW, 

e ^ ^ (2.2) 

x(0) = x~ and x(l) = x^ . 

In this scaling we are studying transitions on a unit time-interval, in which the 
systematic motion of the molecules is large and the thermal noise is of order 1. 
The probability measure vr governing the stochastic boundary value problem 



( |2.2[ ) has density with respect to the Brownian bridge measure ttq arising in the 
case V = 0. The density relating the two measures is found from the Girsanov 
formula, together with an integration by parts (use of Ito formula) and use of the 
boundary conditions on the path 1,15. 16 1, and is given by 



^(x)ocexp(-^j[| G{x;e)d?l (2.3) 



where the path potential G is (with | • | denoting the Euclidean norm) 

e) = ]^\DV{x)\^ - eAF(x). (2.4) 
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Functions are infinite dimensional and tliere is no Lebesgue measure in infinite 
dimensions. Noneth eless , it is instructive to think heuristically of the Brownian 
bridge obtained from ( |2.2j ) with y = as having a probability density (with respect 
to Lebesgue measure) ol the form ifTTl 



exp 







dx 


ds\ 




f 






Jo 


ds 



(2.5) 



(This is the formal limit obtained from the probability density function for a dis- 
cretized Brownian bridge). Let H]_ ((0, 1)) denote the subset of ((0, 1); M^) 
comprised of functions satisfying x{0) = x~ and = x"*". Then, combining 
\23) and \2.5\ , we may think of the probability density for vr as being proportional 
to exp(-2^/e(x) 
is defined by 



^-I^{x)) where the Onsager-Machlup functional 1^ : -?/j_((0, 1)) 



Ie{x) :-- 



dx 



ds 



+ -G{x;e) 



ds. 



(2.6) 



This intuitive definition of Ie{x), via the logarithm of the pathspace probability 
density function, suggests that minimizers of 1^ are related to paths of maximal 
proability. This can be formulated precisely as follows |I3]|4]|. Let B^{z) denote a 
ball of radius 5 in C([0, 1]; M^) centred on z£H\ ((0, 1)) . Then, for any zi, Z2 £ 

^4((o,i)), 

7t(B'{zi)) / 1 



lim 



exp 



Thus, for small ball radius 6, the logarithm of the ratio of the probabilities of the 
two balls is equal to the difference in evaluated at the ball centres. For this 
reason we are interested in minimizers of Jg. 

The following notation will be useful. We define 



Jeix) 

J{x) 



dx 



ds 
dx 



ds 



+ -\DV{x) 



ds, 



ds. 



(2.7) 



Then 



Ie{x) = Je 



AV{x{s))ds. 



(2.8) 



Our goal in what follows is to demonstrate that the F— limit of Is is finite only 
when evaluated on the set of BV functions supported on the critical points of V. 
And, furthermore, that on this set the value of the F— limit is determined by (suit- 
ably rescaled) minima of J, subject to end-point conditions, together with the inte- 
gral of AV{x). Thus we conclude this section with some observations concerning 
minima of J. 
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The Euler-Lagrange equations for the functional J are 

(fx 



ds2 



D^V{x)DV{x) = 0, 



(2.9) 



where D^V denotes the Hessian of V. These equations are Hamiltonian and con- 
serve the energy 

2 



E 



1 



dx 



ds 



1 



\DV{x)\\ 



Thus heteroclinic orbits connecting critical points of V via the equation ( |2.9| ), for 
which E is necessarily zero, satisfy 



dx 



ds 



\DV{x)\. 



(2.10) 



Hence, if there is a heteroclinic orbit connecting critical points of V under either 
of the forward or backward gradient flows 



dx 
d^ 



±DV{x) 



(2.11) 



then this will also determine a heteroclinic orbit in the Hamiltonian system. How- 
ever, the converse is not necessarily the case: there are heteroclinic orbits in the 
Hamiltomnian flow which are not heteroclinic orbits in the gradient dynamics. 

We say that a potential F : M is admissible if F G C3(M^, R) and 

1 . the set of critical points 

£ = {xeM.^ \ DV{x) = 0} 

is finite; 

2. the Hessian D'^V{x) has no zero eigenvalues for every x ^ £; 

3. the weak coercivity condition 



> such that inf \DV{x)\ > 

|a:|>H 



(2.12) 



is satisfied. 



Admissibility implies, in particular, that all critical points for the gradient flows 
( |2.11| ) are hyperbolic. 

For each pair x", x+ ^ £ the set of transition paths is defined as 



X(x-,x+) = {y£ BV{ 



lim y{t) = x^ and y G L^(M) 

t— >±oo 
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Lemma 2.1. For any x ^ X{x , x"^) 

J{x) > \V{x+) -V{x-)\. 



(2.13) 



Furthermore, if the infimum of J over X[x , x^) is attained at x* & X (x , x^), 
then X* G C^(M) satisfies the Euler-Lagrange equation (2.9 1 and has the property 



J{x* 



/oo 
\DV{x*{s))fds. 
-oo 



(2.14) 



Assume that the potential is admissible. If either x or x^ is a local minimum or 
maximum and there exi sts a heterclinic orbi t x* c onnecting x^ and x^ under the 
Hamiltonian dynamics \2.9\ , then x* solves ( |2.1 1[ ) and 



j{x*) = \V[x+) - V{2 



(2.15) 



In par ticular x* is a minimizer of J which attains the lower bound the lower bound 
Proof. Clearly 



J{x) > 



dx 



ds 



> 



> 



— oo 
oo ^ 



— oo 
oo ^ 



DV{x 
{V{x)) ds 



ds 
ds 



ds 



Integrating and using the end point conditions gives the first result. Standard 
regularity results i mply that minim izers x of J are (7^ (M) and s atisfy the Euler- 
Lagrange equation ( 2.9 1. Equation ( |2.14 1 then follows from ( |2.10| ). 

To prove the final result we first note that the linearization of equation ( 2^ i at 
a critical point Xc € £ is given by 



ds2 



D''V{x,fy = 0, 



(2.16) 



whilst the forward gradient flow, given by ( |2.11| ) with a minus sign, has lineariza- 
tion 

dy 



ds 



+ D'V{xc)y = 0. 



(2.17) 



If V is admissible, and Xc is a critical point of V, then (0, Xc) is hyperbolic for 
the Hamiltonian flow and, from ( 2.16 l, has stable and unstable manifolds both of 
dimension N. 

We assume next that x+ is a local minimum of V. Let Ai be the A^-dimensional 
manifold defined by 



M = {(x, -DV{x)) I X G M^} C 



t2N 
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Differentiating equation ( |2.11| ) with respect to t shows that A4 is invariant under 
the flow of equation ( |2.9| ). Since D'^V{x'^) is positive definite there exists e > 
such that all solutions y{t) of equation ( |2 9| ) with the properties |y(0) — x"*"! < e 
and y(0)) = —DV{y{0)) converge to x+ as t tends to infinity. Define next the 
stable manifold 



lim 



}. 



where y solves ( |2.9| ) and satisfies the initial condition y(0) = x, y(0) = v. Since 
the dimension of A4{x^) is the sets Ai{x^) n {|x — < e} and A4n {\x — 
< e} coincide. Thanks to the invariance of A4 it follows that A4{x^) C A4. 
As a result one obtains that 



Jix* 



DV{x*{t))x* = V{x-) - > 0, 



where the last inequality holds because ( 2.11[ ) implies that F is a decreasing Lya- 
punov function along the trajectory x*. Thus, the claim has been proven under the 
assumtion that is a local minimum. The remaining three cases (x+ is a local 
maximum, x^ is a local minimum/maximum) can be dealt with in an analogous 
fashion. □ 

The result shows that only saddle- sadd le connections can give rise to hetero- 
clinic orbits in the Ha miltonian system ( |2.9| ) which cannot be found by in one of the 
gradient flows ( 2. 11 1. As a cl osin g remark in this section we observe that the exis- 
tence of heteroclmic orbit in ( |2.9| ) is a generic property of non-degenerate systems. 
For the Hamiltonian flow the effective dimension of the space is 2 — 1 , because 
it is constrained to a level set of the Hamiltonian. Since the stable and unstable 



manifolds for the Hamiltonian system are both of dimension (cf eq. ( |2.16[ )) we 
expect that, generically, there will be a 1 dimensional manifold connecting any pair 
of equilibria under the Hamiltonian flow. 



3 The r-Limit 



We now determine the F— limit for the Onsager-Machlup function ( 2 6| ). A thor- 
ough introduction to F-convergence can be found in |71 and |6|. ProDlems closely 
related to ours are studied using F-convergence in [ilOllUl . 

We start with some basic definitions and propositions which serve to e xplain 
the form of the F— limit. The key result is contained in Proposition 3.2 which 
demonstrates that every transit ion i s achieved via a finite collection of intermediate 
transitions. The limit Theorem 3.4 is then stated and proved, building on a number 
of lemmas which follow it. 

Recall the definition of admissible potential V, as well as that of a transition 
path. 



Proposition 3.1. Let V be admissible and define for each pair x , x^ 
fiinction $(x~,x"'") by 



£ £ the 



<I>(x" 



inf {J(y) \ yeX{x ,x+)} 
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Then there exists c > such that <I>(x ,x^) > c for all pairs € £ with 

^ x^ . 

We can establish a direct representation of the transition energy $ which avoids 
the usage of infima. 

Proposition 3.2. Let V be admissible and x^ ^ £ be two critical points of V. 
Then there exists a finite sequence {a;i}^=o ^ ^ ^^'^^ ^^^^ = x~, Xk = x+ and 



$(x = ^min{J(?/) \ y £ X{xi^i, Xi)} , 

A straight-forward refinement of the analysis shows that the sequence Xi is 
injective. 

Definition 3.3. Let X be a Banach-space, e > a parameter and : X — )• M a 
family offunctionals. The functional Iq : X — t- M U {+cx)} is the T— limit of as 
e — 7- if for all x £ X and all sequences G X which converge weak- * to x as e 
tends to the liminf-inequality 

liminf > /o(2;) 
e-5>0 

holds, and for all z €z X, there exists a recovery sequence Zs G X which converges 
weak- * to z and satisfies 

limsup/e(z£) < Io{z). 

Theorem 3.4. Let V be an admissible potential. Then the T-limit of the functional 
le as e tends to is 



h{x) 



' Y._H^-{r),x+{T))- j^^V{x{s))ds ^/xeW([0,l]) 

and x{s) G £ a.e. s G [0, 1], 



+00 else, 



where V{x) is the set of discontinuity points of x and x^{t) are the left and right- 
sided limits ofx at t. 



Proposition |3.2| shows that the infimum <I> can be written as a finite sum of 
mini ma. These minima are obtained from evaluation of ( 2.14 i, where x* solves 
( 2.9 ) subject to x*{t) — )• Xi-i/xi as t ^ ±00. Furthermore, if either — 1 or Xi 
is a lo cal ex tremum of V, then x* is also a h eteroc linic orbit in one of th e gradient 



flows ( |2.11[ ), then the minimum is given by d2.15 1. Thus Theorem 3.4 shows that 
the T— limit of can be computed through knowledge of the criticaTpoints of V 
and the set of (Hamiltonian or gra dient ) heteroclinic orbits connecting them. 

We start the proof of Theorem 3.4 by proving a lemma which delivers a lower 
bound for the amount of energy needed to reach one of the stationary points x £ £ 
when starting nearby. We will often write x^ as a shorthand for {x" , x'^}. 
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Lemma 3.5. Let V be an admissible potential. There exists numbers C,eo > 
such that for and all e G [0, £q], ^ £ and all all paths y G X{x^ ,x^) with 
dist(y(0), x^) > e the estimate 

J(y) > 

holds. 

Proof. We first show that there exist numbers 6,C > such that, for all x, y G 
with X ^ £ and with the property — x| < 5, the estimates 

^|y-x|2<|W(y)|2<C|y-x|2 (3.1) 

Indeed, thanks to the smoothness of V and nondegeneracy condition dei{D'^V {x)) / 
for X G <5, for sufficiently small e and large C 

\DV{y) - D^V{x) ■{y-x)\<C\y- xj^, 
\z\<C\D^V{x)-z\ VzGM^ 



hold. Thus 



\DViy)\ > \D'V{x) ■ {y - x)| - C\y - x\' 

> - x\ - C\y - x|^ > - x| 



if |y — x| < min{(5, ^^}- This shows that the first inequality in (3.1 1 holds. The 
second inequality is obtained in a similar way. 

Let y G X{x~ , x"^) and recall that by Sobolev's imbedding theorem y is con- 
tinuous. For each e G (O, min |5, 2^}) we define 

t~^{e) = inm{s €R \\y{t') -x+\ <e foralH'>s}, (3.2) 

and analogous. Then we obtain 

rCfO 1 poo 

J{V) > / \DV{y{s))\ \y{s)\ ds > — (x+ - y) • ds 

Jt+(e) Jt+{e) 

The claim follows since we can derive the same estimate with instead of □ 



Proof of Proposition 3.1 Let r > be the separation of the stationary points of V, 
i.e 

r = min{|x — x'\ | x, x' G <S, x / x'}, 
and let yi G X(x^,x+) be such that limi_^oo J{yi) = ^*(x^,x+). There ex- 



ists a sequence G M such that that dist{yi{ti),£) > r/2 and thus Lemma 3.5 



applied to the translated sequence yi{- — ti) implies that inf/ J{yi) > c, with 
c = 2^min{r/2,eo}^. □ 
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The proof of Proposition 3.2 showing the existence of minimizing connecting 
orbits, is established with the oirect method of the calculus of variations, with the 
aid of the following lemma. 

Lemma 3.6. Letx^ G £. If there exists a minimizing sequence {yi}i^fi G X(x^,x^ 
such that the density pi{t) = 5 + ^^sht in L^(M), then 

there exits a minimizer £ X{x~ , x^) such that J{y^) = x"*"). 

Proof. First we show that the boundary conditions are not lost during the passage 
to the limit. The tightness of pi implies that for each e > there exists (e) G R 
with T~(e) < T+(e) and the property that, for all /, 

1 



lim sup 



{\m\^ + \DV{m{s))f) ds 



< 



AC 



(3.3) 



where C is the constant in Lemma [331 Let 



mm 



{s £R \ \yi{t') - x+\ < £ forallt'>s}, 



and t, analogously. Without loss of generality we assume that t (e) < < tf{e) 



and \yi{0) 



x^\ > 



lim sup 

1^00 2 



e. Then Lemma [33] implies that 
1 



Together with inequality (3.3 1 we obtain that limsup;_^oo */^(^) ^ '^^{^)- Another 
application of this argument with replaced by t~ delivers the desired result 
concerning the boundary conditions: for each e > there exists —00 < t^{s) < 
r+ (e) < 00 such that 



lim sup sup \yi{s) 

Z— >00 S>T + 

lim sup sup \yi{s) 

1^00 S<T~ 



X 



< e. 



(3.4) 



From the bound on iji in (r (e) , r+ (e)) we deduce that the length of the vector 
yi{^) G is bounded. Since J{yi) is bounded the sequence iji is bounded in 
L?{W). We extract a subsequence yi (not relabeled) such that iji converges weakly 
to r/ G and yz(0) converges to z G M^. 

We define the limiting path y^ by 



Equation (3.4 1 implies that lims_j>-|-oo y^{s) 
are satisfied. 



y'^[t) = z + / 77(5) ds 

x^, i.e. the boundary conditions 
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Next we demonstrate that y is a minimizer, i.e. 

lim J{yi) = J(y±). 

t— >oo 

Note that J is weakly lower semicontinuous: the first term is weakly lower semi- 
continuous because it is convex whilst the second term is weakly continuous. To 
see the weak continuity note that 



\DV{y^{s))\ds= lim / \DV{y^{s))\ds 

/R i-oo 
\DV{yi{s))\ds= / \DV{yi{s))\ds. 
-R. >oo J^QQ 



The last equation is again due to the tightness of J{yi), the second but last equation 
holds because of Sobolev's embedding theorem. The weak lower semiconinuity of 
J, coupled with the fact that {yi} is a minimizing sequence, shows that J{y^) = 
lim;^oo J{yi)- 

Finally we show that y^ G L^(M), which entails the claim y^ G X{x^ , 
Define the functional 

and the function g{t) = Gt{y^). Note that (y^(s))s>( minimizes Gt{y) for every 
t, subject to the boundary condition y{t) = y^{t). Testing Gt{-) with a function 
which is affin e on [t, t + 1] and assumes the value for all s > t + 1 together 
with estimate (13. lb delivers the bound 



9{t) < l\y^{t)-x^ 



I' (1 + Cs') = ^1 + f ) \yHt) - (3.5) 



Recall that limi_5.oo y^{t) = and choo se tg > such that \y{t) — x~^\ < Eq 
for all t > Iq, where eq is defined in Lemma 3.5 The function g satisfies for all 
t > to the differential inequality 

— it) > l\DV{y^m' > ^\y^{t) - > ^,g{t), 



with fi = (3^'^)(g - The last inquality is (3.5 I, the penultimate inequlity is due 
to(3J_l. Gronwall's inequality implies the exponential decay bound g'(t) < EoC"^^*" 
Since |y='=(t)p < — 2^(t) we have shown that y € L-'^([to, oo)). The same argu- 
ment also shows that y^ G oo, to]) and thus y^ G L^(M). □ 



The proof of Proposition 3.2 also relies on Lions' concentration compactness 
lemma which we state here for completeness. 
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Lemma 3.7. I\18V Let {pi)i>i be a sequence in L^(M) satisfying 

Pl > in M and lim / pi{t) dt = A, 
Jr 

where A > is fixed. Then, there exists a subsequence (not relabeled) satisfying 
one of the following three possibilities: 

L (compactness) there exists G M such that pi{- — ti) is tight, i.e. 

rti+R 

Ve > 0, 3R < oo such that / Pi{s) ds > A — e; 



ti-R 



2. (vanishing) 

i-t+R 



rt+H 

lim sup / Pi{s) ds = for all R> 0; 
1^°° teM Jt-R 



3. (splitting) there exists < a < A such that for all e > there exists Iq > 1 
and pj , pf G L> (M) such that for all I > Iq 

\\Pi + pf - Pih^ + \\\p}\\l^ -a| + + a- A| < e, 

lim^^^oo dist(supp(/9[),supp(pf)) = oo. 



Proof of Proposition 3.2 Let G X{x , x~^) be a mini mizing sequence of J and 



define pi{t) as in Lemma 3.6 We now use Lemma 3.7 in the following way: we 
show that vanishing cannot occur, and that splitting can only occur a finite number 
of times. 

To rule out vanishing we define, for real number e satisfying 

e G (^0) 2 ~ ^\ \ X, z G £, x ^ z}^ , 

the starting and ending time of the final transition: 

t+ = sup |t G M I inf \yi{s) - x+| ^ ^| > 
t;- =sup{t <t+ I dist{yi{t),£\{x+}) <e} . 
By construction tj" < and Jensen's inequality implies that 

' , , 1 K /dist(x+,£:\{2;+}) -2e\^ , 



(dist(x+,g\{x+})-2£)^ 

^ 2(t^-*r) ■ 
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Since limsup;_^oo f^^Pii^) ds < oo this implies that lim inf ^^i^oo (i;'^ — ) > 0. 
Furthermore, since V admissible there exists a number c > such that inf^^j^- pi ( 

c and thus 

. _ 1 Z"*^ 

limsup(t^ — ti ) < - limsup / Pi{s) ds < oo. 

This pair of inequalities shows that pi does not vanish. 

If the sequence pi spUts, then there exist sequences ai < hi such that [a;, bi] n 
supp(p[ + pf) = ^ and 

lim {ai — bi) = oo 



pi{t) dt > and liminf / pi{t) dt > 0. 

-oo ^^<^ Jbl 



(3.6) 
(3.7) 



We claim now that there exists a third sequence ti G [ai, bi] and a stationary point 
X ^ £ such that lim;_5.oo yi{ti) = x. This is because 



limsup / pi{t) 

(— >cxD J ai 



dt < OO, 



which together with (3.6 1 and the weak coercivity assumption (2.12i delivers that 

lim mi{dist{yi{t),£) \ t G [ai,bi]} = 0. 

l—^oo 

We define next the sequences yj G X{x^ , x) and yf G X{x, x+) as follows: 

r yi{t) if t<ti, 

yht)={ {ti - t + l)yi{tt) + {t - ti)x ifti<t<ti + l, 
y x \ft>ti + l, 

{X {ft<ti-l, 
{t-ti + l)yi{ti) + {ti -t)x ifti-l<t< tu 
yi{t) ift>ti. 



Clearly yj +yf — x e X{x ,x 
that 



and the convergence limj_s.oo yiiti) = x implies 
lim J{yl + yf - x) = lim {j{yl) + J{yf)) = lim J{yi). (3.8) 



Moreover, equation < \3J} implies 



liminf J(y/) > and liminf J{yf) > 0. 

l—>-oo l—^oo 



(3.9) 
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We will show now that x {a;^}- Indeed, if x = x , then yf ^ X{x , x+) and 
by Q 

liminf J{y'i) < lim J{yi) — liminf J(y/). 

Equation ( [3.9[ ) delivers a contradiction to the assumption that yi is a minimizing 
sequence. An analogous argument can be constructed for the case x = x+ and 
thus X € <S \ {x^*"}. 

The construction yields the following equation 

^'(x~, x"*") > lim J(y/ + yf — x) > liminf J{yi) + liminf J{yi) 
>$(x~,x) + <I>(x,x+). 

Progressing inductively we obtain a sequence xo,xi, . . . ,Xfc € £ such that 

xq = x^, Xfc = x+, Xi-i / Xj for alH = 1 . . . A; and 

k 



$(X ,X+) > ^$(Xi_l,Xi). 



Proposition 3.1 implies that k < <I>(x^, x+)/c, placing a finite bound on the num- 
ber of splittmgs possible. The minimizing sequences associated with X{xi-i, Xi) 
are eventually tight and so achieve their contribution to the infimum at a minimizer 
of J. □ 

Proof of Theorem Now the representation formula for the F-limit will be ver- 
ified. Let X G BV{[Q, 1],M'^) be a limit path and x^ G i^^O, 1) such that x^ 
converges to x weak-* in BV{[Q, 1]). First we show that the Laplacian can be 
treated separately. Since weak-* convergence entails boundedness in BV there 
is a constant C > such that limsup^^g l|a;e||L°° < C. Furthermore, thanks to 
Helly's theorem we can select a subsequence which converges pointwise for almost 
every s G [0,1] and therefore by dominated convergence 

lim ! AF(xe(s))ds= ! Ay(x(s))ds. 
<^^o Jo Jo 

Note that only continuity of AV is required for this step. This implies that 

liminf ^(xe) = liminf Je(xe) - /" AV{x{s))ds (3.10) 

£^0 £— !>0 Jq 



with Je as defined in ( |2.7| ). 

Next we consider the case where the set of times s where the limit function 
x(s) G £ does not have full measure, so that /o(x) = oo. Define 

M = meas {{t G [0, 1] | x{t) £}) . 
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There exists 6 > such that 



liminf meas({t G [0, 1] | dist(x(t),^) > S}) > M/2. 



£>0 



Since V satisfies the weak coercivity condition the infimum of exceeds 
some A > 0, on this set with measure exceeding M/2. Thus 



Finally we consider the case where Io(x) is finite. The proof of the liminf- 
inequality follows directly from the definition, as we now show. Since x G BV{[0, 1]) 
and Io{x) is finite, the set of discontinuity points T> C [0, 1] is finite. We assume 
next that V C (0, 1), i.e. x is continuous at the end-points. If x jumps at one or 
both end points the same argument can be repeated with obvious modifications of 
the cut-oils. 

Thanks to this assumption for each r G P the left-sided and right-sided limits 
x'^{t) exist. Since converges weak-* in BV{[0, 1]) to the piecewise constant 
function x Helly's theorem implies that lim.^-^o Xe{t) = x{t) for almost every 
t G [0, 1] \ V. Choose 



and for each r G 2? a pair of such points t G (0, 1) such that t — 6 <t <r< 
<T + S and lime_).o Xe(t^) = x^{t). Setting ye{t) = Xe{et) one obtains 



To justify the boundary condition we modify the path segments ye{t — t / e) in a 
way such that the energy changes only slightly. Define now the function 



liminf Jp(xp) > liminf — 




(5 = I min (r - r' I r, r' G P U {0, 1}, r < r'} . 



> 




x~{t) ifs<^-l 

- s)x-{t) + {s-'i + l)x,{t-) if ^ - 1< s < ^ 

ye{s) if^<S<^, 

[s - ^)x+(r) + + 1 - s)x,it^) if *± < s < ^ + 1 



x+{t) ifs>^ + l. 



The construction above implies that £ X{x {t),x~^{t)) and thus 



J{yl)>^x-{T),x-{T)). 
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Furthermore, the modification of the function affects the energy only in a negli- 
gible way, i.e. 

Hence, the liminf-inequality 

rev 

holds. The same argument also shows that 

liminf Jsixe) = +00, 
e-s>0 

if the set V is infinite. 

To construct a recovery sequence we fix 5 > 0. By definition for each r G I? 
there exists 

yjGX(x-(r),x+(r)) 

such that the support of is compact for all r S P and J(yJ) < <I>(x^ (r), x+(r))+ 
5. It can be checked that for fixed 5 > {) the sequence 

converges weak-* in BV{[0, 1]) to as e — )• 0. 

Moreover, there exists a function d{e) such that lime_>o 5(e) = o(l) and the 
supports of the functions — T)/e), r G P are disjoint as e tends to 0. Hence, 

JeiVsMs)) = E '^-(?^?((- - ^)/^)) ^ E 'J>(^'(t), a:+(T)) + #V6{e) 

rev t£V 

rev 

as s tends to 0. □ 



4 Numerical Experiments 

The aim of these numerical experiments is to illustrate that the F— limit derived in 
sectionJS] accurately captures the behaviour of the problem of minimizing given 
by ( |2.6f when e is small. The intuitive picture of the F— limit is that it is comprised 
of mmimizers with the following properties: 
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1 . The minimize rs ar e B V functions supported on the set of critical points of V 



(see Theorem 3.4 1. 



2. The contribution to the limit functional from the jumps in these B V functions 

(ie from J^) can be expressed in terms of sums of integrals f^°° \DV{x*{s))\'^ds 
where x* is a heteroclinic orbit for the Hamiltonian equations ( |2.9| ) (see 
Lemma [ZT| equation ( |2.14[ )). 

3. Unless both equilibria are saddle points, the n thes e heteroclinic orbits 
will also be heteroclinic for the gradient system ( |2.1 1| ) and then the contri- 
bution to the limit functional from c an b e expressed in terms of sums 



\V{x*{oo)) - y(x*(-oo))| (see Lemma [Zl] equation pl3| )). 



4. Minimizers x* will seek to organize the support of the limiting BV function 
so as to minimize the value of AV{x*){s)ds; this is the second contribu- 
tion to the infimum defined in Theorem 13 .41 

A variety of numerical computations, all of which exhibit these phenomena on 
a range of problems, including high dimensional systems arising in vacancy diffu- 
sion and the Lennard- Jones 38 cluster, may be found in the paper [5 |. The purpose 
of this section is to illustrate the four points above on a single low dimensional 
example and relate the results in an explicit way to the theory developed in earlier 

sections. We employ the potential 1/ : — M given by 

V{XUX2) = {xj + xl) ((Xi - 1)2 + xl) {xl + (X2 " 1)') 

which is shown in Figure [T] The potential has three wells of equal depth, situ- 
ated at Mo = (0, 0), Ml = (1, 0), and M2 = (0, 1). Saddle points exist at Si = 
((2 + \/2)/6, (2 - \/2)/6) and (by symmetry) at 52 = ((2 - \/2)/6, (2 + ^/2)/Q). 
The potential is zero at the minima and attains a value of 2/27 at the saddles. The 
Laplacian of V has value zero at the saddle points, 4 at the minimum Mq and 8 at 
the minima Mi and M2. 

In the followin g nu meric al ex periments we use gradient descent to minimize 
Ie or Je given by ( |2.8[ ) and ( |2.7[ ). In all the experiments we employ a value of 

e^^ = 10~^ which proved to be small enough to exhibit the behaviour of the 
F— limit. We solve the parabolic PDE arising from the L? gradient flow for 1^ 
(resp. Je) by means of a linearly implicit method with stepsize chosen to ensure 
decrease of 1^ (resp. Jg) at each time-step. 

Figure |2] shows heteroclinic orbits for the gradient flow ( |2.1 1[ ) (in green) and 
for the Hamiltonian flow (IS) (in blue). Figure [3] shows minimizers of J^, both 



connecting Mi and M2. It is instructive to compare this figure with the preceding 
Figure[2] The green curve in Figure|3]connects K'li to M2 via Mq and is comprised 
of 4 segments each of whic h has been verified to be approximately given by the 



gradient heteroclinic orbits ( 2.11[ ). The blue curve connects Mi to Af2 via and 



5*2 and has been verifi ed to be approximately given by Hamiltonian heteroclinic 



orbits satisfying \1.9\ . Furthermo re, in the gradient case we have verified that 
the minimizers obey the sum rule ( |2. 15| ) This illustrates the connection between 
minimizers of and solutions of the Luler Lagrange equations for minimizers of 
J, and points 2. and 3. in particular. 

Figures [4] and [5] show minimizers of Jg and 1^ which connect the two saddle 
points Si ana S2 - Several approximate minimizers are shown in each case, found 
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Figure 1 : Contour plot of the potential. The dashed black line is the equipotential 
contour that connects the saddle points, 5i and S2. The three equal-depth minima 
are labeled: Mq, Mi, and M2. 
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Figure 2: The dashed black line is the equipotential contour that connects the sad- 
dle points, Si and 82- The green line is the path the follows the gradient of the 
potential from Mi to Mq via Si and then to M2 via 52- The blue line is a Hamil- 
tonian minimizer (not a gradient flow) that connects Si to 5*2 . 
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Figure 3: Solutions that minimize at e = 10~^. The dashed black line is the 
equipotential contour that connects the saddle points, and 82- The green line 
is the path that starts at Mi and proceeds to M2 via Mq. The blue line is the path 
path that starts at Mi and proceeds to M2 avoiding Mq. 
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from different starting points for tiie gradient flow. The experiments illustrate point 
1. as they show that the solutions concentrate on critical points ofV: in this case 
simply the two saddles. They also illustrate point 4. as they show that, in this case, 
the mmimizers of and are indistinguishable; this is because the Laplacian oi^ 
V is zero at the saddle points. 

Figures [6] and [7| also show minimizers of Jg and which connect the two 
saddle points ana S'2. However the starting points for the gradient flow differ 
from those used to generate Figures ^ and pf in particular they are based on a 
function which passes through the mimmum Mq. As a consequence the minimizers 
also pass through Mq. For there are then multiple approximate minimizers, all 
supported on Si, S2 and Mq. However the support can be organized more or less 
arbitrarily (provided only two transitions occur) to obtain approximately the same 
value of J^; we show a solution where the support is organized symmetrically. 
The situation for Is is quite different: the effect of the Laplacian of V, which is 
4 at Mq and at Si and S'2, means that minimizers place most of their support at 
Mo . The experiments thus again illustrate point 1 . as they show that the solutions 
concentrate on critical points of V. They also illustrate point 4. 



0.0 0.2 0.4 0.6 0.8 1.0 




Figure 4: The xi component of the paths that minimize at e = 10~^. The paths 
start at Si and are conditioned to end at S'2. The ending path depends on the nature 
of the starting path. The results of using three different starting paths are displayed. 
None of the initial paths approach the origin and thus avoid passing through Mq. 
The black curve corresponds to a symmetric minimizer. 

Figures [8] and |9] again show minimizers of and I^, now connecting the two 
minima Mi and M2, and constructed to pass through the other minimum Mq and 
the two saddle points Si and S'2. For there are then multiple approximate min- 
imizers, all supported on the five critical points, one of which is shown in Figure 
[8] a solution where the support is organized symmetrically. The situation for is 
again very different: the eirect of the Laplacian of V, which is 8 at Mi and M2, 
means that minimizers place most of their support at these two points, as shown in 
Figure[9| The single interface in fact contains several transitions, and hence several 
contributions to the F— limit. Furthermore this single interface can be placed ar- 



21 



0.0 0.2 0.4 0.6 O.a 1.€ 




- 0.6 



- 0.4 



- 0.2 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 



Figure 5: The xi component of the paths that minimize at e = 10"'^. The 
paths start at Si and are conditioned to end at 82- The ending path depends on 
the nature of the starting path. The results of using three different starting paths 
are displayed. None of the initial paths approach the origin and thus avoid passing 
through Mq. The black curve corresponds to a symmetric minimizer. Note that 
this plot is indistinguishable from the previous one. 
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Figure 6: The xi component of the path that minimize at e = 10^^. The path 
starts at and is conditioned to end at 5*2. The initial path contains the origin Mq 
and thus the ending path spends an arbitrary fraction of the time at Mq. The ending 
path depends on the nature of the starting path, in particular where the path crosses 
the origin. Only the symmetric minimizer is displayed here. 
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Figure 7: The xi component of the path that minimize at e = 10^^. The path 
starts at Si and is conditioned to end at 82- The initial path contains the origin Mq 
and thus the ending path is dominated by Mq. 



bitrarily; we have shown a symmetric case. The experiments once again illustrate 
points 1. and 4. 
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